Structure of 19 C from Coulomb dissociation studies 
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We investigate the structure of the neutron rich nucleus 19 C through studies of its breakup in the 
Coulomb field of target nuclei. The breakup amplitude is calculated within an adiabatic treatment of 
the projectile excitation, which allows the use of the realistic wave functions for the relative motion 
between the fragments in the ground state of the projectile. The angular distribution of the center 
of mass of 19 C, longitudinal momentum distribution of 18 C and relative energy spectrum of the 
fragments (neutron - 18 C) following the breakup of 19 C on heavy targets at beam energies below 100 
MeV/nucleon have been computed using different configurations for the ground state wave function 
of 19 C. In all the cases, the data seem to favor a 18 C(0 + )®lsi/2 configuration for the ground state 
of 19 C, with the one-neutron separation energy of 0.53 MeV. 

PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.Mn, 27.20. +n 



The neutron rich nucleus 19 C (the last bound odd- 
neutron isotope of carbon) is a strong candidate for hav- 
ing a one- neutron halo structure pLSj, due to its very 
small one-neutron separation energy (S n ) (of the order 
of a few hundred keV). Several measurements reported 
recently do seem to provide evidence in favor of such a 
possibility. The measured longitudinal momentum dis- 
tributions of 18 C fragment emitted in the breakup re- 
action of 19 C on Be and Ta targets at the beam en- 
ergy of 88 MeV/nucleon, and on a carbon target at ~ 
1 GeV/nucleon, have widths (full width at half maxi- 
mum (FWHM)j of 42 ± 4 MeV/c, 41 ± 3 and 69 ± 3 
MeV/c respectively (l]||. The neutron momentum dis- 
tribution, measured in a core- breakup reaction of 19 C 
at 30 MeV/nucleon on a Ta target at GAN1L, shows a 
FWHM of 64 ± 17 MeV/c @. This is about three times 
smaller than that predicted by the Goldhaber model [Q , 
which provides a good description of this quantity in the 
case of stable nuclei. The measured interaction cross sec- 
tion H of 19 C on a 12 C target at beam energy of ~ 960 
MeV/nucleon (1231 ± 28 mb) is seen to be enhanced as 
compared to that of 18 C (1104 ± 15 mb), which signals a 
larger matter radius of 19 C B. Very recent measurement 
of the relative energy spectrum of 19 C in its Coulomb dis- 
sociation at 67 MeV/nucleon on a Pb target |7j shows a 
strong peak at very low relative energy, which is supposed 
to be a characteristic of the halo structure B . 

However, there are still some open issues which make 
the existence of a one- neutron halo structure in 19 C some- 
what unsettled. The detailed structure of this nucleus re- 
mains uncertain, partly because of the large uncertainty 
in its mass Pj9|-p^|. The latest Nubase evaluation ]l3| ] 
gives the value of S n — 160 ± 110 keV. However, analy- 
sis of the data on the interaction cross section of 19 C ||] 
and relative energy spectrum of fragments observed in 
the Coulomb breakup of 19 C JtJ seem to favor & S n — 0-5 
MeV (with a particular configuration for the 19 C ground 
state). Moreover, in a recent measurement of the neu- 
tron angular distribution in the elastic breakup of 19 C 



(in which target nucleus remains in the ground state), 
a broad FWHM of 120 ± 18 MeV/c has been reported 
Ji4| ]. This value is about two times larger than the re- 
sults reported in Ref. wherein a narrow width was 
observed. 

In this paper, we analyze the available data on the 
breakup of 19 C on heavy target nuclei within the frame- 
work of a theory of the breakup reactions which allows 
the use of realistic wave functions to describe the projec- 
tile ground state. By comparing our calculations (per- 
formed with various configurations for the 19 C ground 
state) with the data, we hope to put a constraint on the 
ground state structure of this nucleus. This is expected 
to clarify the issue regarding the existence of a halo struc- 
ture in 19 C since the spin-parity and the configuration of 
the valence neutron play an important role in the forma- 
tion of the halo. We shall consider here only the Coulomb 
breakup process, which dominates the breakup of the 
loosely bound projectiles on heavy target nuclei jl5| in 
kinematical regime below the grazing angle. 

The theory of the Coulomb breakup (CB) of the pro- 
jectile used by us has been described extensively in 
p^pl . The triple differential cross section for the re- 
action, in which a projectile a breaks up into a charged 
core c and a neutral valence particle n on target t, is 
given by 
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Here v a is the a-t relative velocity in the entrance channel 
and p(E c ,Q, c ,Q, n ) the phase space factor appropriate to 
the three-body final state. The amplitude (3fP is given 
by 



tif = (Qn\V cn \^(x ( -\k c );^ n \x (+ \k a )) , (0.2) 
where q n — k n — "fk a , with 7 = m n /(m c + m n ). k n and 
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k a are the asymptotic momenta of neutron and projec- 
tile and m„ and m c are the masses of neutron and the 
core. The first term in Eq. (2) contains the structure in- 
formation about the projectile through the ground state 
wave function <E>^(r), and it is known as the vertex func- 
tion p6[ , while the second term is associated only with 
the dynamics of the reaction, which can be expressed in 
terms of the bremsstrahlung integral |nj . For the expla- 
nation of other quantities in the above equation, we refer 
to@. 

The CB theory is fully quantum mechanical and is also 
non-perturbative. The method retains finite-range effects 
associated with the interaction between the breakup frag- 
ments and includes the initial and final state Coulomb 
interactions to all orders. It allows the use of wave func- 
tions of any relative orbital angular momentum for the 
motion between c and n in the ground state of a. It 
should, however, be mentioned that to obtain Eq. (2) 
fl8|| , it has been assumed that the dominant projectile 
breakup configurations excited are in the low-energy con- 
tinuum (the adiabatic approximation). Furthermore, this 
theory is not applicable to those cases where the valence 
particle is charged. 

Expression for the Coulomb breakup amplitude in the 
factorised form, which also uses the bremsstrahlung in- 
tegral, has been obtained previously fl7f within the dis- 
torted wave Born approximation (DWBA), by making 
the approximation of replacing the vector describing the 
separation of c+n cm. with respect to the target by that 
of c with respect to target in the projectile wave function. 
The important difference between Eq. (2) and the cor- 
responding DWBA expression so obtained is that in the 
latter the vertex function is evaluated at momentum k n 
instead of q n . It has been shown 18 that DWBA with 
this approximation underestimates the {d,pn) breakup 
cross sections at beam energies > 140 MeV. 

In the calculations of the vertex function, we have con- 
sidered the following configurations for the valence neu- 
tron in the 19 C ground state: (a) a lsi/2 state bound to 
a 0+ 18 C core by 0.24 MeV, (b) a ls 1/2 state bound to 
a 2+ 18 C core by 1.86 MeV, (c) a 0d 5/2 state bound to 
a 0+ 18 C core by 0.24 MeV, and (d) a ls 1/2 state bound 
to a 0+ 18 C core with 0.53 MeV. The authors of Ref. § 
use option (d) with a spectroscopic factor (SF) of 0.67. 
The binding potentials in all cases are taken to be of 
Woods-Saxon type having central and spin-orbit terms 
with the radius and diffuseness parameters being 1.236 
fm and 0.62 fm respectively. The strength of the spin- 
orbit term was held fixed to 7 MeV [Q . The depths of the 
central potentials were searched so as to reproduce the 
respective binding energies. The rms sizes of 19 C with 
options (a)-(d) were found to be 3.50 fm, 3.03 fm, 3.00 
fm and 3.23 fm respectively, while the corresponding rms 
sizes of the valence neutron were 9.26 fm, 4.93 fm, 4.43 fm 
and 7.07 fm respectively. The rms size used for the 18 C 
core is 2.9 fm |p0[ . These different wave functions of 19 C 
give rise to different vertex functions and consequently, 



different Coulomb breakup cross sections. 

In Fig. 1, we present a comparison of our CB model 
calculations (performed with configurations (a)-(d) as 
mentioned above) and the experimental data for the an- 
gular distribution of the center of mass (cm.) of the n 
+ 18 C system in the breakup of 19 C on a Pb target at 
the beam energy of 67 MeV/nucleon |Q. The integration 
over the relative energy is performed in the range of - 
0.5 MeV (the same as is done in Ref. 0). We see that 
out of the four cases, the calculation done with option 
(d) (with a SF = 1.0) (solid line) is in the best agree- 
ment with the data. Still, the quality of the fit to the 
data is not as good as that seen in Ref. [7j. This is due 
to the fact that in Ref. |7j] the calculated cross sections 
have been folded with the experimental angular resolu- 
tion. After performing a similar folding, the agreement 
between the CB model results (shown by the dash-dotted 
line) and the data is of the similar nature as that seen in 
Ref. 0. It may, however, be noted that the semiclassical 
Coulomb excitation (SCE) calculations reported by these 
authors |^l| use option (d) with a SF of 0.67, instead of 
1.0 as used by us. With our choice of SF, the SCE theory 
overpredicts the experimental cross sections. 
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FIG. 1. Calculated angular distribution of the centre of 
mass of 19 C in its Coulomb breakup on a Pb target at 67 
MeV/nucleon. The dotted, short-dashed and long-dashed 
curves have been multiplied by 0.2, 100 and 100 respectively. 
The dash-dotted line is the result of the calculation performed 
with option (d) folded with experimental angular resolution. 
The experimental data in taken from Q . 

In Fig. 2, we compare the results of our calculations with 
the experimental data for the relative energy (E re i) spec- 
trum of the fragments for the same reaction as in Fig. 
1. As in Ref. 0, the angular integrations for the 19 C 
cm. are done up to the grazing angle of 2.5°. We note 
that in this case too, the best agreement with the data 
(near the peak position) is obtained with the configura- 
tion (d) (with SF = 1.0) for the 19 C ground state. Here 
again, the SCE calculations done with the same configu- 
ration and SF overestimate these cross sections. A SF of 
0.67 is required to make the SCE results (shown by the 
dash-dotted line) agree with the data in the peak region. 
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One reason for this difference in the CB and SCE results 
could be the fact that in latter the energy distribution of 
the dipole excitation (B(E1)) has been obtained with a 
plane-wave to describe the relative motion between the 
fragements in the final state. Consideration of the nu- 
clear interaction for this may lead to a reduction in the 
cross section as is observed for the case of 8 B |22|. It 
may be noted that in the CB model the final state con- 
sists only of the product of the wave functions describing 
the c+target and n+target relative motions; the wave 
function for the relative motion between the fragments 
does not enter here. We use a plane wave to describe the 
neutron-target relative motion which is valid for the case 
of pure Coulomb breakup of projectiles with a chargeless 
fragment in the outgoing channel. 

We would like to emphasize that, the SCE and CB 
theories use entirely different mechanisms to describe the 
Coulomb breakup process. The CB theory uses the post- 
form scattering amplitude which includes breakup con- 
tributions from entire continuum corresponding to all 
Coulomb multipoles and relative orbital momenta be- 
tween the valence and core fragments. In contrast to this, 
the SCE calculations include contributions from one (as 
is the case in Ref. M) or at the most two multipolarities. 
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FIG. 2. Calculated relative energy spectra in the Coulomb 
dissociation of 19 C on a Pb target at 67 MeV/nucleon. The 
dotted, short dashed and long dashed curves are results of 
multiplication by 0.16, 10 and 10 respectively of the actual 
calculations. The dash-dotted line is the result of the semi- 
classical calculation performed with option (d) (see text). For 
the semiclassical calculation, we use SF = 0.67. The experi- 
mental data have been taken from M. 

Furthermore, often the continuum is curtailed to some 
maximum value. Therefore, the differences seen in the 
SCE and CB models should not be surprising. 

Of course, both the CB and SCE theories underesti- 
mate the relative energy spectrum for larger values of 
E re i. Although, the SCE results are comparatively bet- 
ter in this regard. Improper consideration of the nuclear 
breakup effects could be one of the reasons for this dis- 
agreement. Indeed, Dasso et al. 23 have shown that 



nuclear breakup cross sections dominate the relative en- 
ergy spectrum of fragments for E re i > 0.6 MeV in case of 
the breakup of 11 Be (also a one-neutron halo nucleus) on 
a Pb target at a similar beam energy (72 MeV/nucleon). 
The authors of Ref. § do correct their data for the nu- 
clear breakup effects by scaling the corresponding cross 
sections for the breakup of 19 C on a carbon target. How- 
ever, this is unlikely to be accurate as the scaling pro- 
cedure may not be valid due to the long range of the 
nuclear interaction in the halo nuclei |24| ]. 

The CB and SCE calculations (with configuration (d) 
and SF = 1.0) predict a total Coulomb breakup cross sec- 
tion of 0.78 b and 1.53 b respectively. The latter is even 
larger than the experimental value of the total breakup 
cross section of 1.34 ± 0.12 b (?]]. This again underlines 
the necessity of using a SF of 0.67 in the SCE calcula- 
tions. As far as the CB model results are concerned, one 
should keep in mind that there are non-negligible nu- 
clear breakup cross sections, which must be considered 
together with it for any comparison with the experimen- 
tal total breakup cross sections. In a quantum mechan- 
ical theory, however, one has to add coherently the am- 
plitudes of the Coulomb and nuclear breakup processes. 
Therefore, it is not correct to simply add the cross section 
of the two processes calculated separately for comparison 
with the data on total breakup. Unfortunately, the ex- 
tension of the CB theory to include the nuclear breakup 
effects J2j| is non-trivial and has not been attempted so 
far. 
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FIG. 3. Calculated parallel momentum distributions of 18 C 
from Coulomb breakup of 19 C on Ta at 88 MeV/nucleon. 
The calculations have been shifted to the data to compare 
the widths and the peaks normalized to the peak of the data. 
The experimental data are taken from jy]. 

The parallel momentum distribution (PMD) of the heavy 
charged fragment (which provides an almost unambigu- 
ous information on the existence of a halo structure in 
the projectile 0Jl^,|2^ ^9|) is expected to have a large 
contribution from the Coulomb breakup process in the 
peak region |l|,|l6|. In Fig. 3, we show the comparison 
of our calculations with data ffl| for the PMD of the 18 C 
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fragment emitted in the breakup of 19 C on a Ta target 
at 88 MeV/nucleon beam energy. The results obtained 
with options (a), (b) and (c) are similar to those shown in 
Ref . |30| . The MSU data do not have the absolute mag- 
nitudes of the PMDs. Therefore, we have normalized the 
peaks of the calculated PMDs to that of the data (this 
also involves the shifting of the position of the maxima 
of calculations to those of the data), so that the widths 
of the calculated and the measured distributions can be 
compared. 

The FWHM of the PMDs calculated with the options 
(a), (b), (c) and (d) are 27, 71, 83 and 41 MeV/c re- 
spectively, while that of the experimental one is 41 ± 3 
MeV/c. Thus the option (d) is favoured by this data as 
well. It should be noted that in a study of the breakup 
of 19 C on a Be target J3lj within a core-plus-neutron 
coupling model with a deformed Woods-Saxon potential 
for the neutron-core interaction, an agreement with the 
FWHM of the experimental PMD has been obtained with 
the configuration 18 C(0 + )®ls 1/2 (with S n = 0.5 MeV), 
for the ground state of 19 C. Although this result is for a 
lighter target (Be) , for which Coulomb breakup is not the 
dominant reaction mechanism, yet its comparison with 
that of ours may not be out of the place as the PMD is 
least affected by the reaction mechanism f2?],|3^] . 

In conclusion, we have studied Coulomb breakup of the 
neutron rich exotic nucleus 1 C and compared our cal- 
culations with the existing experimental data in order to 
probe its ground state structure. The breakup amplitude 
is calculated within an approximate quantum mechani- 
cal theoretical model, which assumes that the important 
excitations of the projectile are to the low-energy con- 
tinuum so that they can be treated adiabatically. The 
method permits a finite-range treatment of the projec- 
tile vertex and includes initial and final state Coulomb 
interactions to all orders. 

We find that the calculations performed with the con- 
figuration 18 C(0 + )(g)ls 1 /2 for the ground state of 19 C 
(with a one-neutron separation energy of 0.53 MeV) agree 
with the data on the angular distribution of the 19 C cm., 
parallel momentum distribution of the breakup fragment 
18 C and the relative energy spectrum of 19 C, better than 
those done with other configurations. Therefore, these 
data seem to support this configuration as the dominant 
component in the ground state wave function of 19 C. This 
also gives credence to the existence of a neutron halo 
structure in this nucleus. 

Calculations done within a semiclassical Coulomb exci- 
tation theory (with the same configuration) overestimate 
the magnitudes of the angular distribution of the 19 C 
cm. and the relative energy spectrum of the breakup 
fragments in the peak region; a spectroscopic factor of 
0.67 is required to explain the data within this model. 
The use of a plane wave to describe the relative motion 
between the fragments in the final state in obtaining the 
energy distribution of the dipole excitation used in these 
calculations may be one of the reasons for this overesti- 



mation. It would be worthwhile to redo these calculations 
by considering the nuclear distortion effects in the final 
channel. Proper consideration of the nuclear breakup ef- 
fects in both (semiclassical Coulomb excitation as well 
as Coulomb breakup) models is also necessary to explain 
the relative energy spectrum at higher relative energies. 

The authors would like to thank Takashi Nakamura 
for several helpful correspondences and for performing 
the folding of the calculated cross sections shown in Fig. 
1. 
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